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We present the first observation of the baryon decay A° — > followed by — > pK~iv + in 

106 pb -1 pp collisions at y/s — 1.96 TeV in the CDF experiment. In order to reduce systematic error, 
the measured rate for A° decay is normalized to the kinematically similar meson decay B° -> D + 
followed by D + — > 7r + _K" _ 7r + . We report the ratio of production cross sections (<r) times the 
ratio of branching fractions (B) for the momentum region integrated above pr > 6 GeV/c and 
pseudorapidity range \r,\ < 1.3: a(pp -> A° b X)/a{pp -> B°X) x B(A° -> A+tt-)/S(B -> D+tt") = 
0.82 ± 0.08 (stat) ±0.11 (syst) ±0.22 (B(A+ -> pA"^7r + )). 

PACS numbers: 14.20.Mr, 13.30.Eg, 14.65.Fy 



Weak decays of baryons containing b quarks are a good 
laboratoryfor testing the Heavy Quark Effective Theory 
(HQET) "Pl . The A£ baryon is the ground state of the 
udb quark system, and, in the heavy quark limit, the 
light degrees of freedom are in the state of zero total 
angular momentum Q. Fully hadronic b — > cud tran- 
sitions are more complicated in baryons than in mesons 
because there are diagrams which are not present in the 
decays of the latter. Various extensions of HQET have 



been used to evaluate the A° — > A+7r decay rate 0, 
but the predictions vary over a large range. However, 
in Soft Collinear Effective Theory (SCET) 0, all tree- 
level amplitudes can be properly evaluated, resulting in 
an explicit prediction for the ratio of branching fractions 
B(A° b -> A+tt-)/B(B° -» D+7T") sa 1.7 0. The de- 
cays of A° are also interesting because they may provide 
a means to determine CKM matrix elements with dif- 
ferent systematic uncertainties than the determinations 
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from the decays of B mesons 0. 

This is the first reconstruction of a hadronic decay 
of a b baryon at a hadron collider that does not use 
a J/ip in the final state. In addition, our sample has 
more than an order of magnitude more events than any 
previous sample of fully reconstructed A° decays, and, 
for the same luminosity, is about five times larger than 
a sample of A& — > J/tpA decays. Since b baryons are 
not produced at the B factories operating at the T(4S) 
resonance, studying them comprises a unique facet of 
the B physics program at Collider Detector at Fermi- 
lab (CDF II) (3- In particular, a large sample of fully 
reconstructed A° decays would allow CDF to study other 
properties of b baryons, e.g., to measure the lifetime of 
A°, and also to search for decays of heavier b baryons 
such as — > A^7r ± . 

This paper presents a measurement of a ratio of A^ 
and B° branching fractions multiplied by the ratio of 
production cross-sections, 

*(pp^A° b X) g(A°^A+0 
a{pp^B°X) B{B° -> D+7T-)' {> 

where the a quantities are the cross-sections for A° and 
B° production in the pseudorapidity range |?7| < 1.3 
with momentum in the transverse plane, px, above 
6 GeV/c @. 

The quantity R compares the branching fractions 
of the topologically similar, fully reconstructed decays 
A" — > A+7r~ and B° — » D + ir~ , where the charmed 
hadrons decay via similar three-body channels A+ — > 
pK~ir + and D + — > -k + K~tt + 9]. To extract a branching 
ratio, a good understanding of the absolute reconstruc- 
tion efficiency is needed. However, to obtain R only the 
ratio of the reconstruction efficiencies of two topologically 
and kinematically similar decay modes is evaluated, sig- 
nificantly reducing the systematic errors on the measured 
quantity. 

The upgraded CDF II detector is well-suited for the 
detailed study of weak decays of heavy baryons. In 
particular, the advent of the Silicon Vertex Trigger 
(SVT) |1C|. which uses precise position measurements to 
select events containing weakly decaying heavy hadrons, 
allows CDF II to collect many hadronic decay modes of 
heavy baryons for the first time. This measurement is 
performed using a 106 pb _1 sample of pp collisions col- 
lected by CDF II between February 2002 and June 2003. 
This data sample corresponds to ~ 10 9 b hadron decays 
produced in the central detector region. A full descrip- 
tion of the CDF II detector can be found elsewhere 0- 
The detector components pertinent to this analysis are 
the silicon microstrip vertex detector (SVX II) [II]], the 
drift-chamber central tracker (COT) fl2l | and a three- 
tiered trigger system (Levels 1, 2, and 3). The five 
double-sided layers of the SVX II used in this analysis 
provide up to 10 position measurements. Of these, up to 
five are in the r-(f> Q plane (each precise to about 15 /Ltm), 
three are longitudinal, and two are small-angle-stereo. 
The strips are parallel to the z-axis, longitudinal strips 



are inclined at 90°, and the small-angle-stereo strips are 
inclined at 1.2°. The SVX II spans the radii between 2.5 
and 10.6 cm, and covers the pseudorapidity range \rj\ < 2. 
The COT has 96 measurement layers between the radii 
of 40 and 137 cm. These are organized into alternating 
axial and small-angle-stereo (2°) super-layers. The COT 
has a smaller pseudorapidity coverage (|?7 < 1.3) than 
the SVX II. Both tracking detectors are immersed in a 
1.41 T magnetic field parallel to the z axis. 

This analysis in particular relies on the SVT, which 
operates as a part of the Level 2 trigger system. The 
trigger makes it possible to select events at a rate of 
~ 100 Hz from the ~ 1 MHz interaction rate. The 
components of the three level trigger system pertinent to 
this measurement are the Extremely Fast Tracker (XFT) 
at Level 1 and the SVT at Level 2. The XFT uses 
four axial super-layers of the COT to find tracks with 
p T > 1.5 GeV/c. The SVT combines the XFT measure- 
ment with r — <j) hits from the SVX II detector. The 
track finding is performed using a large look-up table of 
hit patterns. The found track candidates are fitted for 
curvature, angle projected onto the transverse plane and 
impact-parameter |l3j . The impact-parameter measure- 
ment allows the selection of long-lived particles in the 
trigger decision. 

The signal (A°) and normalization (B°) events are col- 
lected using the same trigger. At Level 1, two tracks 
must satisfy pt > 2.0 GeV/c, a scalar sum of trans- 
verse momenta pri + PT2 > 5.5 GeV/c, and an an- 
gular separation projected onto the transverse plane of 
Ac6 < 135°. At Level 2 the transverse momentum cuts 
are repeated, and it is required that each track has im- 
pact parameter do > 120 /im, with an angular separation 
between the tracks projected onto the transverse plane of 
(2° < Atfi < 90°). Finally, the distance evaluated in the 
transverse plane from the primary vertex to the two-track 
intersection point must be greater than 200 /im. 

Additional criteria are imposed on the triggered sam- 
ple in order to reject as many background events as pos- 
sible while keeping most of the signal. To reconstruct 
A° — > A+7r~, minimum acceptable COT and SVX hit 
requirements are imposed, and all valid four track com- 
binations are considered. The px of the proton candi- 
date from the A+ and the ir~ candidate from the A° 
must be greater than 2.0 GeV/c, which strongly favors 
these particles to be the two which caused the event to 
pass the trigger. The pt of the proton candidate must 
be larger than the pr of the ir + candidate from the A+ . 
The pr of the A° and A+ candidates must be greater 
than 7.5 GeV/c and 4.5 GeV/c, respectively. In Equa- 
tion (Q), cr A o and ergo are defined for pr > 6 GeV/c. The 
events in the data sample must satisfy pr > 7.5 GeV/c, 
and the difference is accounted for by using the Monte 
Carlo simulation based on the pt distributions of both 
A" and B° measured in data [l4| . 

Each of the unstable particles (A°, B°, A+, and D + ) 
is reconstructed by considering all valid combinations of 
tracks and requiring them to satisfy the decay hypothe- 
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sis. The charmed hadrons A+ and D + are reconstructed 
first: each triplet of tracks that satisfies the selection cri- 
teria (detailed below) is constrained to pass through the 
same point, called the decay vertex. The decay vertex is 
determined by varying the track parameters of the stable 
daughters within their uncertainties to minimize the x 2 . 
The A+(Z) + ) candidate is then combined with a fourth 
track to form a A®(B°) candidate. The full topology of 
the decay is then imposed in another kinematic fit, re- 
sulting in a simultaneous measurement of the A° (B°) 
and A+ (D + ) vertices. 

Using these measurements, the reconstructed invari- 
ant mass of the A+ must be between 2.269 GeV/c 2 and 
2.301 GeV/c 2 . Other selection criteria rely on L xy , the 
projection onto the x-y (transverse) plane of the decay 
length measured from the production vertex to the decay 
point; the production vertex is estimated by the position 
of the beamline averaged over each run calculated for the 
z coordinate of the secondary vertex. A product of the 
proper decay time and the speed of light, ct, is also used. 
It is derived from L xy : ct — (L xy ■ Pt)[tiic/pt), where 
L xy is the decay vector of A° (B°) projected onto the x- 
y plane, pt is the transverse momentum, pt is the unit 
vector in the direction of the transverse momentum, and 
m is the world average mass of the A° (B°). In order 
to suppress the combinatorial background from the in- 
teraction point, we impose ci(A°) > 225 /im (compared 
to the b baryon mean decay length of 368 ± 24 /an [l5|l. 
Calculated relative to the A° decay point, ct(A+) must 
be ci(A+) > —65 ^m. For a true A+, a small negative 
ct may arise due to resolution effects; on the other hand, 
the ct(A+) of combinatorial background candidates may 
have large negative values. The distance of closest ap- 
proach in the transverse plane of the trajectory of the A° 
candidate to the primary vertex must be less than 85 fim. 
The x xy °f the A° and A+ kinematic fits are required to 
be less than 30 and 20 respectively, where x xy 1S a X 2 ~like 
goodness-of-fit quantity using only the track parameters 
in the transverse plane. 

The normalization mode (B° — » D + tt~) is recon- 
structed using selection criteria identical to those of the 
signal mode, except for a different invariant mass re- 
quirement for the D + candidate and no analogy to the 
Pt(p) > Pt{^ + ) cut. The D + candidate invariant mass 
must be between 1.848 GeV/c 2 and 1.888 GeV/c 2 . The 
distributions of the invariant mass of A^ — > A+7r~ and 
B a — > D + tt~ candidates are shown in Figures ^ and 
respectively. 

Figure shows the binned likelihood fit to the invari- 
ant mass distribution of A° candidates. The large Gaus- 
sian peak at 5.6 GeV/c 2 is the signal. The dash-dotted 
curve corresponds to the exponential combinatorial back- 
ground. This component is constrained by the data in 
the invariant mass region above the A° mass. The small 
asymmetric peak at 5.5 GeV/c 2 (solid line) corresponds 
to contributions from fully reconstructed i?-meson de- 
cays resulting in a final state with four tracks, where at 
least one track is misidentified. This shape is obtained 



Other B-meson decays 
■ Other A b decays 

Combinatorial background 
- Four particle B-meson decays 




A+tt Mass (GeV/c ) 

FIG. 1: Aj — > A+7r~ yield with binned likelihood mass fit. 
The background shapes are defined in the text. 



■ Physics backgrounds 

■ Combinatoric background 




DV Mass (GeV/c ) 

FIG. 2: 5° -* D+tt~ yield with binned likelihood fit. The 
background shapes are defined in the text. 



using a full detector simulation of these modes. It is con- 
sistent with the shape of B° — > D + ir~ candidates found 
in the A° sample. The dotted and dashed curves cor- 
respond to all the other B-meson and A° backgrounds, 
respectively. These shapes are determined from a large 
parametric Monte Carlo sample which includes all known 
decays of B + , B°, and B s and A° hadrons. Finally, there 
is a very small Gaussian distribution (not shown) from 
the Cabibbo-suppressed mode A° — > A+K~ (fixed to an 
expected 8% of the signal yield [l6|)- The total distribu- 
tion is Ftot = G s ignal + -Ecomb + -Ffour —track + -Foth er— B ~r 



F. 



other-A? 



+ G 



ATK- 



where G indicates a Gaussian dis- 



tribution, E indicates an exponential, and F indicates a 
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more complicated functional form. 

In the fit, the width of A b signal is fixed to 
26.4 MeV/c 2 , obtained by scaling the B° width in data 
by the ratio of widths from the Monte Carlo simula- 
tion. The relative contribution of each background type 
in the fit is guided by two constraints: the first describes 
the normalization of -Ffour-track relative to -F thcr-B (i.e., 
Nfour-track/Nother-B), where N is number of events; the 
second describes the normalization of ^other-A relative 

to (Fothcr-B + ^four-track), (is., -^other A? / (-^four-track + 

-^other-B))- The value of each constraint is inferred from 
the relative abundance of the background types in the 
large parametric Monte Carlo sample. The value of 
the (iVf our _t r ack/-^Vother-B) constraint is checked by recon- 
structing the B° —> D + n~ mode among A+7r~ candi- 
dates from the region of the invariant mass correspond- 
ing to the Ag signal. The total x 2 of this fit is 80.6 for 88 
degrees of freedom, corresponding to the fit probability 
of 70%. 

Figure [3 shows the mass fit for the D + ir~~ candidates. 
The large Gaussian peak at 5.27 GeV/c 2 is the B° signal. 
The dashed curve corresponds to the exponential combi- 
natorial background. The dotted curve corresponds to 
the background from decays of other b hadrons. The 
total x 2 °f this fit is 70.9 for 94 degrees of freedom, cor- 
responding to the fit probability of 96%. 

The quantity R defined in equation Q is calculated 
from the signal yield according to 

N B0 B(A+^pK-7r+) e(A°)' [) 

where the first factor is the ratio of observed signal 
yields, the second factor is the ratio of the daughter 
branching fractions [l^ , and the third factor is the ratio 
of reconstruction efficiencies calculated from the Monte 
Carlo. The signal yields are 7V A o = 214 ± 19 and 
Nga = 790 ± 32 respectively. Each reconstruction ef- 
ficiency is defined for the A"(_B°) with pr > 6 GeV/c 
and |tj| < 1.3. The exact configuration of the CDF II 
detector varied over the course of collecting the data 
used in the analysis. However, the ratio of reconstruc- 
tion efficiencies is stable within statistical errors across 
the different periods of running, with the average value 
e(B°)/e(Ag) = 1.65 ± 0.03. From Equation © we thus 
obtain R = 0.82 ± 0.08(stat). 

The systematic uncertainty on the measurement of R 
is dominated by the error on S(A+ — > pK~n + ), yielding 
a relative error of 27% 0. Since this uncertainty is in- 
dependent of our measurement, it is quoted separately. 
Other sources of systematic uncertainty are the imper- 
fect knowledge of the A° lifetime (1^%), the production 
p T spectra of both A° b (7.6%) and B° (4%), the A° b po- 
larization (7%), and the A+ resonant substructure (1%). 



The uncertainty due to the finite size of the Monte Carlo 
samples is 1.9%. The uncertainty due to the difference 
between the proton and tt + trigger efficiency is 0.6%. The 
systematic uncertainties due to the background shapes in 
the invariant mass fits are 4.3% for the A b and 0.9% for 
the B°. These uncertainties arise from the lack of de- 
tailed knowledge of a variety of branching fractions con- 
tributing to background shapes that are obtained from 
the Monte Carlo simulation and fixed in the fit. To eval- 
uate the uncertainty due to these shapes, the branching 
fractions of the largest decay modes contributing to each 
of the shapes were varied simultaneously in the simu- 
lation, and the shapes were reevaluated. Uncertainties 
on the mass resolutions of both A b and B°, which are 
also fixed in the mass fits, contribute 2.8% and 1.8%, re- 
spectively. Finally, the contribution of the A" — > A^K~ 
shape is varied by a factor of 2, contributing 1.6% to the 
systematic error. The total systematic error excluding 
the uncertainty on £>(A+ — > pK~ir + ) is 13.5%. 

A direct comparison with a theoretical prediction of 
Rbr = B(A° b -> A+n-)/B(B° -> dW+tt") « 1.7 i 
can be performed if one assumes that cr(Ab) and cr(B Q ) 
have the same dependence on px, and then use /baryon/ fd 
from high-pT measurements. From Ref. we obtain 
/baryon//d = 0.25 ±0.04, which yields R B r m 3.3 ± 1.2, 
in agreement with j5[. 

In summary, we have observed the decay A® —> A+7r _ 
for the first time, and measured R — 0.82 ± 0.08(stat) ± 
0.11 (syst) ± 0.22 (B A +). The overall error is dominated 
by the large uncertainty on £>(A+ — > pK~ir + ). The A b — > 
A+7r~ sample is the largest 6-baryon sample in existence, 
and, once augmented by new data, can be used for a 
variety of other A^ measurements, including its lifetime 
and production properties. 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nuclcarc; the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the 
Republic of China; the Swiss National Science Founda- 
tion; the A. P. Sloan Foundation; the Bundesministerium 
fur Bildung und Forschung, Germany; the Korean Sci- 
ence and Engineering Foundation and the Korean Re- 
search Foundation; the Particle Physics and Astronomy 
Research Council and the Royal Society, UK; the Russian 
Foundation for Basic Research; the Comision Interminis- 
terial de Ciencia y Tecnologfa, Spain; in part by the Eu- 
ropean Community's Human Potential Programme un- 
der contract HPRN-CT-2002-00292; and the Academy 
of Finland. 



[1] For a review, see A. F. Falk, "Introduction to hadronic B at an Asymmetric B Factory', |arXiv:he p-ph/9812217 

physics," Chapter 2 of 'The BaBar Physics Book: Physics 



7 



and references therein. 

[2] A. F. Falk and M. Neubert, Phys. Rev. D 47, 2982 (1993). 

[3] X. H. Guo, Mod. Phys. Lett. A 13, 2265 (1998). 

[4] C. W. Bauer, S. Fleming and M. E. Luke, Phys. Rev. 
D 63, 014006 (2001); C. W. Bauer, S. Fleming, D. Pir- 
jol and I. W. Stewart, Phys. Rev. D 63, 114020 (2001); 
C. W. Bauer and I. W. Stewart, Phys. Lett. B 516, 134 
(2001); C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. 
Rev. D 65, 054022 (2002). 

[5] A. K. Leibovich, Z. Ligeti, I. W. Stewart and M. B. Wise, 
Phys. Lett. B 586, 337 (2004). 

[6] I. Dunietz, Z. Phys. C 56, 129 (1992). 

[7] D. Acosta et al. Phys. Rev. D 71, 032001 (2005) 

[8] CDF II uses a cylindrical coordinate system with the z- 
axis along the nominal beam-line. The transverse plane 
(a;, y) is perpendicular to the z-axis. Azimuthal angle, 
<j>, is measured from the :r-axis. Polar angle, 8, is mea- 
sured from the z-axis. Pseudorapidity is defined as r\ = 
tanh - (cos#). Transverse momentum, pr, is the com- 
ponent of the particle's momentum projected onto the 



transverse plane. 
[9] Inclusion of the respective charge conjugate modes is as- 
sumed throughout this paper. Specifically, AjJ — > A^ty + , 
B° -> D~tt + , A' -> pK+Tv~, and D~ -> -k-K+tt'. 
[10] W. Ashmanskas et al, Nucl. Instrum. Meth. A 518, 532 
(2004). 

[11] A. Sill et al, Nucl. Instrum. Meth. A 447, 1 (2000). 
[12] A. Affolder et al, Nucl. Instrum. Meth. A 526, 249 
(2004). 

[13] The curvature of a particle track is inversely propor- 
tional to the transverse momentum (pr)- The impact- 
parameter of a particle track is defined as the distance 
of closest approach of the particle track to the primary 
vertex in the transverse plane. 

[14] S. Yu, PhD thesis, University of Pennsylvania, 2005 (un- 
published). arXiv:hep-ex/0504059 . 

[15] S. Eidelman et al, Phys. Lett: B 592, 1 (2004). 

[16] This choice is motivated by K. Abe et al, Phys. Rev. 
Lett. 87, 111801 (2001), however we assign a large sys- 
tematic error to it. 



